VHE 7-ray Emission from Passive Supermassive Black Holes: 

Constraints for NGC 1399 
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ABSTRACT 

Very high energy (VHE; >100 GeV) 7-rays are expected to be emitted from the vicinity of 
super-massive black holes (SMBH), irrespective of their activity state. In the magnetosphere 
of rotating SMBH, efficient acceleration of charged particles can take place through various 
processes. These particles could reach energies up to i? ^ lO^^cV. VHE 7-ray emission from 
these particles is then feasible via Icptonic or hadronic processes. Therefore passive systems, 
where the lack of a strong photon field allows the VHE 7-rays to escape, are expected to be 
detected by Chercnkov telescopes. We present results from recent VHE experiments on the 
passive SMBH in the nearby elliptical galaxy NGC 1399. No 7-ray signal has been found, neither 
by the H.E.S.S. experiment nor in the Fermi data analyzed here. We discuss possible implications 
for the physical characteristics of the system. We conclude that in a scenario where particles are 
accelerated in vacuum gaps in the magnetosphere, only a fraction ^ 0.3 of the gap is available 
for particle acceleration, indicating that the system is unlikely to be able to accelerate protons 
up to £; - IQi^eV. 

Subject headings: gamma rays, galaxies: individual (NGC 1399) 



1. Introduction 

Spheroidal systems (such as elliptical galaxies, 
lenticular galaxies, and early- type spiral galaxies 
with massive bulges) are commonly believed to 
host super-massive black holes with masses in the 
range Afa ^ = (10* ^ - 10^) Mp, in their central re- 
gion (e.g. 



Richstone et al.lll998l : lFerrarese fc Ford 



20051 ). Many SMBH are active and give rise to 
emission over a wide range of frequencies. Some 
have been detected from radio to high energy 
7-rays (- 1 GeV). The SMBH is essential for 
the production and stability of rclativistic jets. 
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In many objects the inner, compact jets emit 
VHE 7-ray emission through various leptonic 
and/or hadronic processes. In blazar-type ob- 
jects, where the jet is pointing towards the ob- 
serve r, VHE 7-ray emission is frequ ently observed 
re.g.. lHinton fc Hofmannll2009l: l^kes 2009,). Of- 
ten, this can be successfully modeled by non- 
thermal electrons upscattering either their own 
synchrotron radiation or upscattering an exter- 
nal photon field, although a hadronic origiii (e.g., 
lAharonianl I2OOOI: iMiicke fc Protherod l200lh can- 
not be excluded. In these cases the gamma-ray 
emission is released through the jets which in turn 
are powered by the SMBH. 

VHE emission may also be indirectly related to 
SMBH independent of any electromagnetic activ- 
ity. The presence of a SMBH steepens the po- 



tential well and hence the density profile of dark 
matter in the central regions of galaxies. The 
rate of annihilation of Dark Matter particles (e.g. 
Berezinskv et al.lll994l) will thus be increased re- 
sulting in enhanced emission of the gamma-rays 
generated in the annihilation. 

Several models for the direct production of 
7-ray emission in the v icinity of SMBH have 
been proposed (see e.g., lAharonian et al. 19841: 



Mastichiadis fc Prothero3 Il990t jSlane fc Wagh 



1990l:lB oldt fc Gos h'1999': Lcvinson2000':'Ncrono v et al 
2004: Neronov fc Aharonian 2007: Ricger fc Aharonian 



20081: llstomin fc Solll2009l : IOsmanovl2010HLevinson fc Rie ' 



20111 ). Possibly, some of these mechanisms could 
also be responsible for acceleration of cosmic-rays 
to energies E ^ 10^^ cV and beyond. In or- 
der to avoid attenuation on circumnuclear fields, 
VHE photons could escape only if the SMBH does 
not produce too much low-energy radiation. The 
SMBH would have to be passive at low energies, 
i.e. most of the radiative losses would have to 
occur at high energies. In all cases a large mass of 
the central object is an important characteristic 
for generating a high VHE fiux. 

Correlations involving the SMBH masses and 
properties of their host galaxies have been in- 
vestigated by many authors. In particular Mbh 
is found to be linked to the central stellar ve- 
locity dispersion (e.g., ICebhardt et al.l I2OOOI ) or 
to the mass of the host galaxy bulge (e.g., 
iMagorrian et al.lll998 ). These observational scal- 
ing laws, in addition to confirming the ubiquity 
of SMBH, have suggested that their activity is 
tightly l inked to the evolution o f their host galax- 
ies (e.g. lFerrarese fc Fordl l2005M. During the early 
stages of galaxy evolution SMBH accrete matter at 
high rates and are observed as bright QSO (Quasi- 
Stellar Objects). Even if such systems generate 
VHE gamma-rays, the dense photon fields asso- 
ciated with the Quasar phase would pair-absorb 
the VHE radiation. The average radiative out- 
put at low photon energies (e.g., in the optical 
band) decays from rcdshift z > 3 to z = by 
almost 2 orders of magnitude. The majority of 
SMBH in the local universe are hosted in sys- 
tems of low accretion rate and are therefore not 
embedded in dense radiation fields. In passive sys- 
tems (i.e., SMBH hosted in nuclei without bright 
signatures of broad-band activity and very low lu- 
minosity in longer wavelengths). If VHE gamma- 



rays, if generated, can escape from the nuclear 
region without suffering from strong absorption 
via photon-photon pair absorption. 

While the detection of VHE gamma-rays in 
blazar-type systems is facilitated by the superlu- 
minal motion (the apparent luminosity is boosted 
and the optical depth r elated to photon-p hoton 
absorption is reduced, see lSchlickeiserl ( 19961 )). this 
is not the case for non-blazar systems. Hence, 
proximity and low luminosity in the IR/optical 
domain increase the possibility of a detection in 
the VHE band. Observations of passive systems 
geac i hence contribute to our understanding of the 
physics and properties of galactic nuclei. In addi- 
tion they might give us an insight on Ultra High 
Energy Cosmic Ray (UHECR; E > Ax lO^* eV) 
sources. 

Here we present GeV limits and discuss obser- 
vations of the passive SMBH in the core of NGC 
1399, the central galaxy of the Fornax cluster, that 
were conducted with the H.E.S.S. Cherenkov tele- 



scope array (jPedaletti et al.ll2007l ) and discuss its 
implications for gap- type particle acceleration and 
emission models. 

2. The test object NGC 1399 

The most massive black holes can be found 
in the large ellipticals at the center of galaxy 
groups and clusters. The nearest SMBH reside 
in the giant radio galaxy Centaurus A and in the 
central region of the Virgo cluster, especially in 
M87. These galactic nuclei have already been de- 
tected in VHE observations by, e.g., the H.E.S.S. 
experiment. The detected VHE emission might 
well be assoc iated with features of their power- 
ful inner jets (JAharonian et al.ll2006l . l2009al) . Also 
SgrA*, in the galactic center, would be an obvious 
good candidate, despite its low black hole mass of 
Mbh ~ 4 X IO^Mq, thanks to its proximity. How- 
ever, it is very close to other established Galactic 
gamma-ray emitting sources. Hence, these galax- 
ies do no represent the best candidates for the 
studies described above. 

The most nearby cluster visible from the south- 
ern hemisphere is th e Fornax cluster. It w as ob- 
served by H.E.S.S. (JPedaletti et al.ll2007l ). The 



giant elliptical galaxy NGC 1399 is located in the 
central region of the Fornax cluster at a distance 
of 20.3 Mpc. A SMBH of Mbh ^ 5.1 x IO^Mq re- 



sides in its central region (jGebhardt et alj 120071) . 
The nucleus of this galaxy is well known for 
its low emissivi t y at all wavelengths, e.g., see 
O'Connell et al.l ( 2005 ) and references therein. 
The galaxy shows low-power antiparallel jets, 
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with a luminosity of 

Hz and 10^° Hz Ixnieen et al 



erg s 



between 10^ 
19881 ). The out- 



eration is (Levinson 2000) 



(1) 



flows are confined in projection within the optical 
extension of the galaxy. The jets are initially tran- 
sonic, then decelerate in the inter-stellar medium 
and end in lobe -like diffuse structu re at ~ 9 kpc 
from the center ( Killeen et al.l 19881 ). X-ray images 
reveal cavities suggesting that the radio-emitting 
plasma is producing shocks in the hot X-ray emit- 
ting gas (IShurkin et al.l 120081) . The estimated jet 
power is Ljet — 10^^ erg s~^. 

Given the suggested low photon density at low 
energies, the high mass of the central SMBH, the 
absence of luminous jets and proximity of the 
galaxy, NGC 1399 emerges as the best candidate 
for a study of passive SMBH at very high energies. 

3. Estimate of VHE luminosity from Pas- 
sive SMBH 

The expected luminosity and energy range of 
VHE 7-ray emission from magnetospheric pro- 
cesses in SMBH can be estimated as follows: 

If the central black hole is spinning and accretes 
matter that carries magnetic flux, it will develop 
a rotating magnetosphere. This magnetosphere 
may be similar to those of neutron stars. Ambi- 
ent charged particles can be accelerated to very 
high energies in the magnetosphere and radiative 
cooling will result in high energy radiation. A full 
treatment of the acceleration and of the resultant 
VHE 7-ray spectra requires detailed modelling of 
charged partic le trajectories and in teractions, as 
shown in, e.g.. lNeronov fc Aharonian (2007) . 

It is however possible to estimate the expected 
VHE 7-ray luminosity based on some simple ar- 
guments: 

The maximum available pot ential due to field line 

rotation is AF ~ (a/2) thB ( Thorne. Price fc Macdonaldl 
19861 ). where a = J/Jmax (with Jmax = GM"^ /c 



If particles are accelerated along the magnetic field 
lines, the primary radiative loss is curvature emis- 
sion (note that photo-meson production is highly 
inefficient in low ambient radiation fields): 



P„ 



2 g'c 4 

KJ J- ^'r■^^r\r 



(2) 



with q the particle charge, -Rcurv the curvature ra- 
dius and 7 the Lorentz factor of the particle. Bal- 
ancing acceleration in the gap by losses implies 



qAVe 



.jh/c. 



(3) 



The expected emission spectrum from curva- 
ture radiation then extends up to VHE energies 
and the energy of the curvature photons does not 
depend on the mass of the emitting particle, i.e.. 



ii;,,_ = ||E2_ ^ 4.8 a^/^My^B?/^ 



hV^URr 



th 



9 ^4 

1/2 



TeV, (4) 



where A/g == Mbh/IO^Mq, B^ = B/W G, and 
fg = 1.5 X lO^^Afg cm is the gravitational ra- 
dius. Hence, for realistic magnetic field strengths 
{B » IG), curvature emission is expected to 
peak in the Fermi and/or H.E.S.S. energy range. 
In the case of curvature-limited losses, an acceler- 
ated proton could reach (provided the potential is 
large enough) ultra-high energies of 



E^ curv ^ 2 X lO^^a^/^Ml'"^ 



1/4 



hR. 
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curv 



rV* ^^urv ^y^ 
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(5) 

In terms of the Goldreich- Julian number den- 



sity ncj = Q,Bcos9/{2'n:cq) (JGoldreich fc Julian 
1969| ). the maximum VHE power of the gap for 

(rie/ncj) < 1 is given by 



1^^ — 



the maximum angular momentum) is the black 
hole spin parameter, rn denotes the radius of the 
event horizon and B the strength of the ordered 
magnetic field component. Provided a gap of size 
h exists, the effective potential for particle accel- 



1] 



TIq J curv ^ * 
f ne\ ngB 



(6) 



(gAK: 
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where fl — a (c/2r/f ) is the angular frequency of 
the black hole and t] ^ 1 is a geometrical factor. 



and where the volume element has been approxi- 
mated by dV ~ irhrjj. Using Eq. ^ one finally 
obtains 



^. - h'(^)^'(-] rlc 



ncj 



8 X m^'^T] a 
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V'^GJ 



Bi Mi 



erg s 



(7) 



Note that for 



riGj, ^7 is roughly a frac- 



tion (h/rg)'^ of the maximum Blandford-Znajek jet 
power Lbz- 

In the case of NGC 1399, the black hole mass is 
Mq ~ 0.5. The spin parameter a depends on the 
formation history of the black hole. If its growth is 
dominated by gas accretion and the accreted mass 
is of the order of the initial mass of the black hole, 
then it is possible to reach spin s approaching a ~ 1 
(see e.g.. lVolonteri et al.ll2007l ). This is believed to 
be true for the SMBH hosted in elliptical galaxies 
in the nearby universe, but might be vio lated if 
chaotic accretion occurs (JKing et al.ll2008l ). 
A lower limit on the magnetic field {B ^ 10 G) 
in the inner disk might be obtained by assuming 
equipartition with the observed radiation fields. 
However, as the radiative output is very low, the 
disk is unlikely to be of the standard-type. Thus, 
a more realistic value is given by equipartition 
with the accretion energy density, i.e., B^/8tt = 
l/2p{ro)v'^{ro), where p is the mass density and 
Vr is the radial infall velocity. The plasma density 
at the Schwarzschild radius is p{rs) ~ 3 x 10~^^ 
g/cm'^, when ex trapolated frorn the value at the 
accretion radius (|Pellegrinill2005l) , assuming a free- 
fall profile. With Vr{rs) = c/-\/2, this would sug- 
gest a field strength B(rs) ^ 1300 G. A similar 
value is obtained if the inner disk would be of the 
ADAF/radiatiyely in efficient accretion flow-type 
(jNaravan et al.lll998l) . In fact, due to its low over- 
all luminosity NGC 1399 has become a prom inent 
ADAF candidate source fe.g.. lNaravanl2002l) . One 
thus expects the accretion rat e in NGC 1399 to 
be close to to ~ 10~'*TOEdd (cf. iLoewenstein et al 



20011 : lNaravanll2002l) and the inner disk magnetic 
field to be close to i34 ~ 0.1. Note that these val- 
ues are compatible with the requirement that the 
jet power (Ljot — 10*^ erg/s) is ultimately pro- 
vided by a Blandford-Znajek-type process. 

The above suggests that the VHE luminosity in 



NGC 1399 could be as high as 



L^ ~ 10«a2 (^ 



V«GJ/ yrg 



erg s ^ , (8) 



and hence high enough to be detectable by a VHE 
telescope array like the H.E.S.S. experiment. 

4. Radiation Fields of NGC 1399 

In constructing the spectral energy distribution 
(SED) of NGC 1399, observations were selected 
in order to consider only its nuclear region. To 
achieve this, at any wavelength only the observa- 
tions taken with the smallest aperture reported in 
the literature were included (see Fig. [T]). In most 
energy bands the smallest apertures exceed the 
projected linear scales considered above by sev- 
eral orders of magnitude. The measurements are 
hence to be considered as upper limits of the fiux 
emitted from the vicinity of the SMBH. 

4.1. Gamma- ray observations 

In the VHE band NGC 1399 has been observed 



with the H.E.S.S. instrument (Pedaletti et al 



2007, 2008) 



The observations did not result in 
a detection of an unresolved source and an upper 
limit (99.9%) on the isotropic VHE 7-ray lumi- 
nosity of 



L^{> 200GeV) < 9.6 x 10^° 



erg s 



(9) 



has been derived (jPedaletti et al.ll2007l ). Extended 
emission would not be expected given that the an- 
gular resolution of H.E.S.S. exceed s the angular 



size o f the host galaxy (~7 arcmin, iJarrett et al 
(J2OO3I ) ) . The given upper limit assumes a spectral 
index of F = —2.6 for the photon spectrum. 

While the H.E.S.S. upper limits refers to the 
integrated flux above 200 GeV, lower-energy 
gamma-ray emission can be studied using data ob- 
tained with the LAT instrument onboard Fermi. 
The data-set consists of 23 months of public Fermi 
data (04 August 2008 - 07 July 2010). The data 
set has been analyzed with the public software re- 
leased by the Fermi Collaboration (ScienceTools- 
v9rl5p2-fssc-20090808). Only events between 200 
MeV and 100 GeV were considered and standard 
selection cuts for point-sources (as outlined in 
the Cicerone Fermi manual) were applied. The 
spectral analysis has been performed through an 



unbinned likelihood fit (gtlikc tool in the pub- 
lic software). First a likelihood fit on the entire 
data range was made with a power law functional 
form dN/dE = NoiE/Ep)-'^, where Nq is the 
normalization at the pivot energy Ep and F is the 
spectral index. The object is not significant on 
the whole energy range (2.4cr). The background 
has been modelled as being composed from the 
sources in the first year Fermi Catalog and from 
two diffuse backgrounds. The diffuse backgrounds 
(both galactic and extragalactic) have been mod- 
elled from the ones supplied along with the public 
software. Upper limits have been derived at a 
95% level for the three energy bins that are of 
most interest for the models investigated in this 
paper. The assumed spectral index for the upper 
limits is F = 2. The bins are l-3GeV, 3-lOGeV, 
10-lOOGeV. Results are shown in Fig. [H 

4.2. Potential Absorption of VHE 7-rays 

The absence of gamma-ray signals might, in 
principle, result from an annihilation of gamma- 
ray emission along the line-of-sight towards the 
observer. In the case of NGC 1399, however, any 
such absorption is unlikley: The cross section a^^ 
of this process depends on the product of the en- 
ergies of the colliding photons. For VHE photons, 
the most effective interaction at the peak of the 
cross section is with background photons of en- 
ergy 

ew l(S/lTcV)"^ eV. 

The optical depth resulting from this absorption, 
in a source of luminosity L and radius R, is given 

by 

L (e) cT^^ 



Energy (eV) 

10"' 10* 10-^ r 10" 10* 10' ro' 10'° 10" lo'lo" 



T (E, R) 



(10) 



AnRec 

For NGC 1399, the visibility of a 200 GeV photon 
would require i (e = 5eV) < 4 x 104i(i?/100 r^,) 
erg/s. This condition is satisfied (Fig. [T]) since 
the flux measured within a small aperture is very 
low. If, instead, more energetic hard photons are 
considered, E > ITeV, the constraints are placed 
by photons in the near- infrared regime (e = leV). 
In the absence of measurements with small aper- 
tures, (weak) constraints can only be derived from 
the fairly high near-infrared flux obtained with 
an aperture which is likely to be dominated by 
starlight. This would set a lower limit on the size 
of the emitting region of i? ~ 3 x lO'^rg in order 
to allow VHE photons to escape. 
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Fig. 1.— The SED of NGC 1399. All archival data 
are for the core region. The archival points are: VLA 
radio data f black: ISadler at al.l (| 19891 '): core component 
<0.4" ); 15/im ISO IR data (red square; iTemi et all 
(|2005l ): 5" aperture, no host subtrac t ion); 2MASS J- 
band data (red star: ISkrutskie et al.l (|2006l l: 4" aper- 
ture, no host subtraction); HST optical data (green; 
lO'Connell et all (|2005h : 0.2" ape rture), and Chan- 
dra X-ray upper limit (solid line; iLoewenstein et al.l 
(|200ll ): 3" aperture). The purple upper limits (95%) 
are the Fermi data presented in this paper. The 
blue line is the H.E.S.S. spectrum 2a upper limit 
(|Pedaletti et al.ll2007l ). The thin black line represents 
a toy curvature spectrum, i.e., a delta-function popu- 
lation with energy as in eq. (|4]) and with maximum lu- 
minosity from eq. ((8} assuming a fully-developed gap. 



An additional photon fleld that might possi- 
bly interact with VHE 7-rays produced near the 
SMBH is pro vided by synchrotron radiation emit- 
ted by pairs (jStawarz fc Kirkll2007[ ). However, the 
critical luminosity for this to occur is well above 
the estimated VHE output, so that this process is 
unlikely to limit the escape of VHE photons. 

Lastly, absorption on the diffuse extragalactic 
background light does not cause any significant 
attenuation due to the proximity of NGC 1399. 

We conclude that photon-photon pair absorp- 
tion is expected to be low throughout the gamma- 
ray range considered here and should not signifi- 
cantly affect any potential VHE emission emitted 
from the vicinity of the SMBH of NGC 1399. 



5. Discussion and Conclusion 

If efficient gap-type particle acceleration and 
curvature emission occurs in NGC 1399, VHE 7- 
ray emission should have been detected by Fermi 
and/or H.E.S.S. According to eq. ([8]), the non- 
detection of NGC 1399 thus either suggests that 
(i) the strength of the ordered, magnetospheric 
field component is only a small fraction (< 0.1) 
of the disk magnetic field value estimated above, 
(ii) on average only a small fraction of the gap 
{[h/rgY ^ 1) is available for particle acceleration 
and/or (iii) the charge density n^ in the vicinity of 
the black hole is much smaller than the Goldreich- 
Julian density, i.e. (ne/"-Gj) *C 1. 

Option (i) would be i ncompatible with th e jet 
power Ljet ^ lO'*^ erg/s (jShurkin et al.ll2008t ) be- 
ing provided by a Blandford-Znajek-type process. 
Option (ii) would require that the charge den- 
sity around the black hole exceeds the Goldreich- 
Julian density ncj '^ 2 x IQ-'^{B/WG) cm^^ 
This could be the case if efficient pair pro- 
duction occurs near to the black hole (e.g.. 



Moscibrodzka et al.l2011tlLevinson fc Riegeii2011l) 
The accretion rate for NGC 1399 inferred above 



(cf. also lNaravanll2002f ) seems indeed close to the 
critical value rhc ~ 2 x 10~^ where annihilation of 
MeV photons in a two-temperature ADAF could 
lead to th e injection of seed charge s with density 
Tie > riGj (jLevinson fc Riegeiil2011l ). If the accre- 
tion rate would be sufficiently high, i.e., m > mc 
ensuring Uf. > ncj , a substantial part of the gap is 
expected to be screened, suggesting {h/vg)'^ <C 1. 
This could make the anticipated VHE output, 
eq. (|5]), consistent with the VHE upper limits de- 
rived above. However, it would also imply that the 
available potential, eq. ([1]), is reduced down to a 
level where UHE proton acceleration will no longer 
be possible. On the other hand, if the accretion 
would be such that rh < rhc, implying ng_ < hqj 
(option [iii]), fully-developed gaps {h ^ Vg) may 
exist. This could again make L^ consistent with 
the observationally inferred VHE upper limits. 
However, an additional plasma source would then 
be needed to establish the force-free outflow be- 
lieved to be present on larger scales. One plausi- 
ble scenario relates to pair cascade formation in 
charge-starved magnetospheric regions due to the 
absor ption of inverse Compton up-scattered pho- 
tons (jLevinson fc Riegeiil2011l ). When applied to 



NGC 1399, the estimated multiplicity appears in- 
deed sufficiently high (Af >, 10^) to allow the pair 
density to approach the Goldrcich- Julian density. 
As the optical depth for pair production across 
the gap is larger than unity even if the gap is not 
fully restored (i.e., for h < Vg), cascade formation 
can occur on scales < rg and thereby limit the 
potential gap size. Detailed modelling would be 
required to self-consistently calculate the gap size 
and associated VHE emission. On the other hand, 
the estimated VHE 7-ray output L^ is roughly 
a fraction (h/rgY of the maximum Blandford- 
Znajek jet power Lbz- Assuming iez ~ -^jot, the 
VHE flux upper limits imply that (JijrgY < 0.1, 
again yielding conditions not conducive for effi- 
cient proton acceleration beyond 4 x 10^^ eV. 

Nearby, passive supermassive black holes are of- 
ten believed to be prime candidates for magneto- 
spheric, gap- type particle acceleration and emis- 
sion scenarios. The non-detection of VHE gamma- 
ray emission in NGC 1399 by current VHE instru- 
ments now suggests that this object is unlikely to 
be an efficient UHECR proton accelerator. 
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